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We investigate the baryonic contributions to the dilepton yield in high energy heavy ion collisions 
within the context of a transport model. The relative contribution of the baryonic and mesonic 
sources are examined. It is observed that most dominant among the baryonic channels is the 
decay of A'^*(1520) and mostly confined in the region below the p peak. In a transport theory 
implementation we find the baryonic contribution to the lepton pair yield to be small. 

I. INTRODUCTION 



One of the main goals of relativistic heavy ion physics is to explore the possibiUty of studying strongly interacting 
matter in a highly exotic form: the quark gluon plasma. In this context the production of dileptons in high energy 
heavy ion collisions has received particular attention in recent years. This interest is due to the fact that dileptons, once 
^ ' produced, essentially decouple from the strongly interacting system and will reach the detectors mostly unscathed. 
' Hence the dileptons serve as a promising probe, as their production rate is biased towards regions of high densities 

and high temperatures. Information on those is only available indirectly if one is limited to hadronic observables. 
' It is possible to roughly separate the dilepton invariant mass region in three parts: the "soft" dileptons essentially 
, lie below the (f) peak and the "hard" dileptons lie beyond the J/V"- The region in-between has been coined the 
"intermediate mass region". It is important to realize that different physics can be at work in those different regimes, 
, mainly reflecting the different epochs of the relativistic nuclear collisions. For example, the hard dileptons will receive 
an important contribution from the Drell-Yan process which will happen in the first instants of the interaction. In 
this work, we shall restrict our attention to the soft part. This region is of interest as it has been the focus of 
detailed recent experimental measurements. Those were done by the HELIOS/3 § and CERES g collaborations, 
O respectively. Those heavy ion experiments have reported an excess of lepton pairs over measurements involving 
^ . proton- nucleus collisions. This has stimulated a great deal of theoretical activity |^-||] which has been summarized 
■ recently §§■ 

The theoretical interpretation of the dilepton excess has mostly concentrated on possible in-medium effects. One 
school of thought attributes the low mass abundance to a dropping vector meson mass, precursor of a chiral symmetry 
restoration in the dense medium [^. Another evaluates the in-medium p spectral function and finds it considerably 
broadened (its peak value is not shifted appreciably), owing mainly to the coupling of the p with baryonic resonances 
1^, with the N*(1520) playing a dominant role as first pointed out in To be precise, the effect is distributed over 
several channclsj IR] . However, the current experimental data are also consistent with a scenario without in-medium 
modifications |11 I4]. Certainly, data with better statistics and an improved mass resolution are needed to extract 
any possible in-medium effects as well as to distinguish between the different in- medium modifications. 

Several issues deserve further investigations in this matter. The dilepton calculations based on in-medium spectral 
functions receive a substantial contribution from the baryonic sector ^ . This feature is at first view puzzling, since the 
baryons are a minor component of the hadron population at CERN energies ]l3[ . A rough estimate of the contribution 
of baryons to the dilepton channel had been provided in , and turned out to be about a factor of two below the 
contribution of the cj-Dalitz decay. With this in mind, we have set out to explicitly investigate the role of baryons in 
dilepton production in relativistic nuclear collisions. Our paper is organized as follows: we first explicitly examine the 
role of the iV*(1520) resonance, which will turn out to be the largest baryon dilepton channel. We also examine a few 
subtleties associated with its off-shell behaviour and its coupling to electromagnetic radiation. We then compare rate 
calculations with each other in order to get a feeling for the magnitude of the different contributions in a somewhat 
idealized thermal environment. We then model the dynamics of the full nuclear collision in the UrQMD approach 
1^^ . The transport model will cover aspects that are related to the possible importance of pre-equilibrium dynamics, 
while insuring a proper reproduction of the hadronic observables. It is important here to insist on the following: 
throughout this work the effective spectral function of the vector mesons are the vacuum ones , i.e. we will work with 



1 



an unmodified pion form-factor. Of course the spectral function associated with the electromagnetic current-current 
correlator will receive contribution from the medium, such as the Dalitz decays of mesons and baryons inside the 
fireball [|l2|. Therefore in this language the low mass lepton pairs come from hadron reactions and Dalitz decays. 
While formally the lepton pair production process can be linked to the imaginary part of the vector meson self-energy 
I psf , its connection to Dalitz decay is clear Bear in mind that a consistent transport treatment of in- medium 

spectral functions is a topic that still requires development. In keeping with our focus on the role of baryons we then 
formulate a prediction for the dilepton spectra measured in recent low-energy runs of the CERN SPS. We finally 
conclude. 



II. BARYONIC INTERACTIONS 



Since one of the main issues we wish to address in this work is the role played by baryons in the production of low 
mass lepton pairs, this section first deals with technical issues that arise in such calculations. The baryons will mostly 
manifest themselves through their radiative decay channel into a dilepton. To fix the ideas, we explicitly consider 
the case of spin 3/2 baryon resonances which will turn out to be the most important contribution in any case. The 
channel we are considering is thus R — > Ne^e^ , where R denotes the baryon resonance. The interaction Lagrangian 
is 

r-RNj = ^v^^^e(zi)^.7ArT3VArF''^ - ^v^ge„^(z2)T3ra.VjvF''^ + /i.c. , (2.1) 

where &fj,,y{z) = g^v — 1/2(1 -I- 22)7^7,^. F is either 1 or 75 depending upon the parity of the resonances and T is a 
3/2 — > 1/2 isospin transition operator or a 2 x 2 Pauli matrix, depending on the isospin of the resonance. and V' 
correspond to Rarita-Schwinger and nucleon spinors respectively, F^^'^ is the electromagnetic field tensor, and M is 
the nucleon mass. The influence of the off-shell parameter z is seen in calculations of electromagnetic transitions of 
nucleon resonances into different multipolarities [ p7[ . It should be mentioned here that when the spin 3 /2 particle is 
on-shell the results are independent of parameter z. Explicit calculation shows that the Rarita-Schwinger projection 
operator 



A,.(g) = {^ + Mn){-g,. + + -7.7. - (2-2) 



which appears in the squared matrix element for the decay of the 3/2 resonance contracted with z-dependent vertex 
terms 7^ vanishes for on-shell particles. In other words, only the first term of the remains operative for vertex 

involving on-shell 3/2 resonances. The presence of the second term in the Lagrangian is important in order to keep 
the electric quadrupole (E2) and magnetic dipole (Ml) transitions independent |L7| . Quantitatively however, we find 
that the contributions from the second term to the dilepton yield are smaller than that from the first one by an order 
of magnitude. 

Before going into more detailed discussions, we point out some differences between this work and some current 
popular approaches that rely on a nonrelativistic (NR) reduction of 

C = Z^vl/^^-TF^,^^ (2.3) 
rrip 

to calculate the p spectral function. To cast it into the form relevant for the NR case, one neglects the lower component 
of the Dirac spinor and upon simplification one obtains 



C = 



'^^'^^ ^RiSkPkL^ - poSkqk)ipN (2.4) 



where Sk {k = 1, 2, 3) is the 3/2 spin transition operator |Q, p^ = {po, pk) is the p meson field, and q^ = (i^Tqk) is 
its four momentum. 

The approach represented by Eq. (|2.lD enables a study of the influence of the off-shell parameter z. This is done 



later on. By setting z = —1/2 and 92 ~ and using the the Vector Meson Dominance model (VMD), Eq.(2J) can be 
cast in the form of the Lagrangian of Eq. ■ We note that it is satisfying to have a context where the importance of 
relativistic effects can be appreciated quantitatively. Those are found to be of some importance in the lower invariant 
mass region, as we will see shortly. 
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FIG. 1. Solid and dashed lines represent results obtained with the non-relativistic and relativistic versions of Eq. (2.3). The 
dashed-dotted line shows the result obtained with Eq. 



Let us consider the A^*(1520) resonance. In the interaction Lagrangian of Eq. ( |2.lD we use parameters obtained in 
Ref. p9[ |, where a fit to pion photoproduction multipoles was performed. To be quantitative, we take gi = —1.839, 
.92 = 0.018, zi = —0.092 and Z2 = —0.024. There exists other parameter sets which yield comparable ^^"^ ''^^ 

the photoproduction data is concerned, but in addition we will require sensible results for the calculation of radiative 



decay channels. The above quoted values of the parameter set yields the following widths: F 



Af*(1520)^Ar7 



0.78 MeV 



and r^.(i52o)^Arp — 22.35 MeV. The experimental values are 0.55±0.1 MeV and 24.0 ± 6.0 MeV respectively. Using 
VMD with the nonrelativistic interaction Lagrangian Eq. (2.4) and fitting its only constant in order to repro duc e 
F^v^Afp, one obtains fp — 5.5 and Fjv^ = 0.88 MeV. Using the same fp in the relativistic Lagrangian, Eq. (^.3|), 
produces Fjv^ = 1.13 MeV and Fjvp = 30 MeV. 

Given a Lagrangian, the differential partial width into a lepton pair of invariant mass M is calculated by standard 
techniques [Q. We first investigate the effects of the different interactions. We fix the mass of the A^*(1520) at its 
central value {i.e. at the peak of the vacuum spectral function). The results are shown in Fig. |^. The solid and 
dash ed curves labeled Non- Rela tivistic and Relativistic I respectively correspond to the cases where the Lagrangian 
(2.3) and its NR reduction (2.4) have been used with the same coupling constant, fp — 5.5. Comparing those two 
curves thus permits a direct assessment of relativistic effects. Those are found to be largest at the photon point and 
they become smaller at larger invariant masses, where the nucleon has a smaller momentum. The relativistic and NR 
treatments should then be in agreement at the larger invariant masses and this is indeed the case. The dashed-dotted 
curve (Relativistic II) represents the differential width obtained using Eq. (2.1). The difference between the two 
relativistic results appear to be a simple overall scaling. We have verified that this overall shift in fact corresponds 
to the ratio of the two radiative decay widths in the two models. We recall the versatility and the success of the 
relativistic approach represented by Lagrangian ( |2.l| ) in terms of its capability to reproduce pion photoproduction 
multipoles and decay widths into and pN final states. However, the quantitative differences with a nonrelativistic 
treatment are overall not large. The relativistic approach can also lend itself to a study of off-shell effects (owing to 
its ^-dependence): a topic towards which we now turn. 

We have calculated dT/dM"^ for iV*(1520) Ne+e' in three different ways. First using the on-shell mass ("On- 
Shell", in Fig. |^), then integrating over the vacuum spectral function of this resonant state with ("Off-Shell I"), and 
without ("Off-Shell 11") the z-dependent interaction. This latter case was achieved by setting z = —1/2. It is clear 
that off-shell effects can arise not only from the mass distribution but also from the associated vertex containing the 
factor &p,i,{z). A comparison of those three approaches is shown in Fig. ||. It is seen that off-shellness can cause 
differences of at most a few percent. Summarizing this section, one concludes that in this case off-shell effects are not 
numerically important and that relativistic effects are slightly larger. However, a nonrelativistic treatment remains a 
sensible approximation. 
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FIG. 2. The solid line shows results when A''*(1520) is on-shell. The dashed-dotted and dotted curve depicts the results 
where the differential width has been integrated over TV* (1520) spectral function with (Off-shell I) and with out off- shell 
interaction {z = —1/2) respectively (Off-shell II). 



Before going into the transport calculation and thereby deaUng with possible compUcations owing to the dynamics, 
we present thermal rate calculations of various lepton pair sources and compare them with with each other. As 
discussed above, the relativistic effects are not large although not necessarily negligible. Note that the difference 
between the relativistic and non-relativistic results might become substantial, if evaluated in a different frame of 
reference |2^. However, in the present work we will employ the calculated decay within a transport model, where the 
widths are always evaluated in the rest frame of the resonance and then simply boosted to the appropriate matter 
f ram e. Thus, in this context using the non-relativistic expression for the widths is a good approximation. Using Eq. 
( ^.4| ) and VMD, the differential partial width for the decay of baryonic resonance Ri into a nucleon and a lepton pair 
of invariant mass M can be written as |lCf| 



III. RATES 




(3.1) 



where 



kp for p 
= 2kl + SM^ 



wave 



for s — wave. 



(3.2) 



and F{kp) — A^/(A^ + fc^), with A = 1.5 GcV. In the above, Jr-np is the coupling coupling constant, St is the spin 
sum, and kp is the momentum of the rho meson in the rest frame of the baryon resonance. 
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FIG. 3. Relative contributions of the baryonic resonances to the dilepton production rate and their sum. 

We should note that the vector dominance couphng used here for the A^*(1720) and A(1905) over-estimates the 
measured photon decay width ^. Thus our results for the dilepton rates due to these resonances are in fact upper 
limits. 

The p mass distribution is 



(3.3) 



One uses the differential decay width to derive the differential dilepton production rate. For a general process 



a b + e+e-, this is ||20(1 

dRa-*b+e+e- 



dM2 



(27r)2 



dAP 



dEaPafaiEa) / dx[l + ^{6^)] 



(3.4) 



where Et = {EaE^ +paPlx)/ma and E^ — (m^ + 1711 ~ M^) / C^rria) j^^]. In the above we have assumed that a and 
b were fermions and thus the /'s are Fermi-Dirac distribution functions. 

We first calculate and compare dilepton rates associated with the radiative decay of several hadronic resonances. 
The species we consider and the associated rates are shown in Fig. ^. We show representative results obtained at 
normal nuclear matter density. Adding to the sum of baryonic contributions the signal from tt — tt annihilation and 
that from uj decay one obtains the net rates shown in Fig. ^. 
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FIG. 4. Comparison of the rates obtained in this work with that of others at T = 150 MeV. The short-dashed curve 
corresponds to tt-tt annihilation only, and the fuU, long dashed and dashed dotted curves refer to the rates calculated here 
(DGKK), by SYZ and by RUBW respectively. See the main text for references. 
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FIG. 5. The relative importance of mesonic and baryonic contribution to the dilepton rate is shown. The sohd hne corresponds 
to the total yield while the dashed line represents the sum over baryonic contributions. Dashed and dotted curves refer to the 
TT-TT annihilation and uj Dalitz decay, respectively. 

The two mesonic sources we have considered are actually the largest ones in the invariant mass range we have 
chosen to consider. Also shown on Fig. || are the rates obtained by Steele, Yamagishi, and Zahed (SYZ) |g2| and by 
Rapp, Urban, Buballa, and Wambach (RUBW) |^^. At very low baryonic densities, the three rates are quite close 
to each other. As the density grows our rates remain quite close to that obtained in the SYZ formalism. At the 
highest baryonic density we have considered, our results approximately lie in between those of SYZ and RUBW, for 
low invariant masses. A marked difference exists in the results of RUBW, that can be understood in terms of the 
spectral function re-summation technique leading to a suppression at the p peak To provide a reference point, the 
analysis of particle ratios for SPS-coUision at 160 GeV per nucleon gives a chemical freeze-out temperature of about 
T = 170 MeV and a baryon density of about p = po (see e.g. j25|] ). Finally, the relative importance of our mesonic and 
baryonic rates is shown in Fig. |^. The baryons dominate the low invariant mass region while the tt — tt channel picks 
up near the p peak. Representative results at normal nuclear matter density are shown. The relative importance of 
mesons and baryons does not change considerably if one uses the chemical freeze-out temperature of 170 MeV. It is 
finally important to realize that the actual contribution from the w Dalitz decay to the final measurements is much 
higher than it appears from Fig. p[ the dominant part is due to the decay of the final state w's. 



IV. TRANSPORT CALCULATION 



The nuclear dynamics and dilepton production is studied here within the framework of the Ultra-Relativistic 
Quantum Molecular Dynamics model, UrQMD [Q. In addition to earlier molecular dynamics models its collision 
term describes the production of all established meson and baryon resonances up to about 2 GeV with all corresponding 
isospin projections and antiparticle states. 

In the present model, light meson formation at low energies is modeled by multi-step processes that proceed via 
intermediate heavy meson and baryon resonances and their subsequent decay |l^ . The pole masses, decay widths 
and branching ratios of all resonances are taken from |^^. However, due to the experimental uncertainties a certain 
range of parameters might be used to obtain a consistent fit to cross section data. As an example, the production of 
UJ mesons is described in the UrQMD model by the formation and the decay of the A^*(1900) resonance. It decays in 
35% of the cases into Ntt and 55% into Nui. In line with data at SIS energies (1 AGeV), the rj production proceeds 
not only via A^*(1535), but invokes also nucleon resonances with masses from 1650 MeV to 2080 MeV. A detailed 
description of the resonance formation cross sections and comparisons to data is give in Ref. |l^. For earlier works 
on dileptons within the same framework see e.g. ||27|] . 
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Dilepton production is calculated in the following way. We consider two distinct classes of processes. First, direct 
decays of vector mesons, i.e. the p and uj mesons. Second, the Dalitz decay of mesons and baryons during the reaction, 
specifically the w and ai meson as well as the iV*(1520). The Dalitz decay of the pion and eta is treated at the end of 
the simulation according to Ref. ||^. By treating the direct decay of the p we implicitly also include the contribution 
form the pion annihilation channel (see e.g. [^). In order to avoid double counting, we do not allow p mesons created 
via the decay of the A^*(1520) or ai-meson to decay into dileptons, since we treat the Dalitz decay of these states 
explicitly. Since in UrQMD the w-meson does not have a p — tt decay channel, there is no double counting in this 
channel. One should note at this point that treating the Dalitz decays as a two-step process via an intermediate p 
meson can lead to erroneous results in a transport description: in this case lepton pairs below the minimum mass of 
m < 2m,r cannot be produced while this is of course kinematically allowed for a Dalitz decay. 

In order to extract the dilepton yield from UrQMD, we extract for each of the aforementioned channels the lifetime 
and four-momentum. The properties of the individual hadrons are given by the collision dynamics. If a meson is 
produced in an elementary collision (may it be due to string fragmentation, meson or baryon resonance decay or 
annihilation) it gets a mass according to a Breit-Wigner distribution and a 4-momentum vector according to the 
kinematics of this single reaction. The meson is then assumed to travel on a straight line trajectory until it (i) decays 
or (ii) collides with another particle and its 4-momentum changes. Any change of the 4-momentum of the meson is 
considered as a destruction of the original meson and the creation of a new meson with new quantum numbers and a 
new 4-momentum. Thus for each hadron under consideration UrQMD provides a list of four-momenta q and lifetimes 
At, {{g, At}i}^, where the index h labels the hadron, and i labels the different four-momenta and time-intervals for 
which the hadrons exists during the course of the collision. The resulting dilepton spectrum is then obtained from 

^?^.i:(^'.r„^,,-(M).(v^-M))^_^^ (4.1) 

for the the direct decay of vector-mesons. For the Dalitz decays, on the other hand we have 

dM ^ V ' 

If absorption is negligible, this approach is equivalent to the method of adding one dilepton (with appropriate 
normalization) at each decay vertex. However, this "shining" method gives better numerical statistics compared to 
a single decay of the vector meson. To get an upper limit on the possible dilepton radiation, particles are allowed to 
couple to the dilepton channel also during their formation time. If we insist that the particles shine only after their 
formation time is over, we have verified that the net rate is only reduced by « 25%. 





M [GeV] 

FIG. 6. Invariant mass spectrum for all pair momenta, for Au-Au collisions at 40 AGeV. Data are from |^ and are shown 
only to set a scale. See the text for details. 
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Since we are specifically concerned with the role of the baryons in the production of lepton pair, we present a set of 
predictions for recent CERN low-energy runs where baryon density effects should be larger than for the high energy 
runs. Those predictions appear in Figs. ||, ^ and ^. 
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FIG. 7. Same caption as for Fig. |6|, but only for low pr lepton pairs. 
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FIG. 8. Same caption as for Fig. but for high pt lepton pairs. 

To provide a visual scale reference, we also show data by the CERES collaboration for 160 AGeV collisions ||]. Aside 
from the clearly visible w-peak, we predict that there will be little difference between the "high-energy" (160 AGeV) 
and "low-energy" (40 AGeV) dilepton invariant mass spectra. At an invariant mass of 400 MeV, the 160 AGeV 
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calculation is 20% higher (including the CERES normahzation) than the 40 AGeV one Here, we have worked 
with a mass resolution of AAf/M = 1%. If this resolution can be achieved in experiment, the omega-peak should 
be clearly visible, providing a strong constraint for presently available model calculations. Furthermore, we find that 
the contribution of the A^*(1520)-Dalitz (thick dashed curve) is very small, even at these low bombarding energy. We 
finally note, that a calculation based on the model of [|j with the addition of the A^*(1520) channel and constrained 
by the final state hadronic spectra as predicted by the RQMD model for a 40 AGeV collision, gives virtually the same 
results. This model also is in fair agreement with the dilepton spectra measured in the higher energy collisions ||T^. 



V. CONCLUSION 



We see that our transport results lead to a very small baryonic contribution to the net dilepton yield for the recently 
done low energy heavy ion runs of the CERN SPS. This is to be contrasted with the baryon influence on the net rate, 
as shown in Fig. ^. The reason for this difference is two-fold. First, the baryonic contribution are severely cut down 
by the CERES acceptance which we apply to our transport results in Figs. ^ and 0. Second, the baryon content of the 
central rapidity region is smaller than its mesonic content, owing to the dynamics of the stopping power systematics 
and of the phase space for particle creation. The small baryon contributions obtained here can be compared with 
similar results obtained with a hydrodynamic model p9| . Therefore it really does appear that the ideal test to pin 
down the relevant dilepton production scenario would involve a low energy run. A high resolution measurement at 
the positions of the vector meson vacuum masses would highlight the high baryon density. In such an environment 
the collective behaviour of the many-body calculations leads to a suppression of the vector meson peak In this 
context we have generated a set of predictions and the corresponding dilepton results of the CERN low energy runs 
are eagerly awaited. We are also looking forward to HADES measurements at the GSI where it has been shown 
that the dilepton contribution from N*(1520) Dalitz decay is most important in the invariant mass region 0.35 < M < 
0.75 GeV/c2 |l|. 
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